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Polydisperse fluids are encountered everywhere in biological and Industrial processes. These fluids 
naturally show a rich phenomenology exhibiting fractionation and shifts in critical point and freezing 
temperatures. Here, we study the impact of size polydispersity on the basic nature of Lennard-Jones 
(LJ) liquids, which represent most molecular liquids without hydrogen bonds, via two- and three- 
dimensional molecular dynamics computer simulations. A single-component liquid constituting 
spherical particles and interacting via the LJ potential is known to exhibit strong correlations 
between virial and potential energy equilibrium fluctuations at constant volume. This correlation 
significantly simplifies the physical description of the liquid, and these liquids are now known as 
Roskilde-simple (RS) liquids. We show that this simple nature of the single-component LJ liquid 
is preserved even for very high polydispersities (above 40% polydispersity for the studied uniform 
distribution). We also investigate isomorphs of moderately polydisperse LJ liquids. Isomorphs are 
curves in the phase diagram of RS liquids along which structure, dynamics, and some thermodynamic 
quantities are invariant in dimensionless units. We And that isomorphs are a good approximation 
even for polydisperse LJ liquids. The theory of isomorphs thus extends readily to multi-component 
systems and can be used to improve even further the understanding of these intriguing systems. 


I. INTRODUCTION 

Polydisperse fluids are ubiquitous in both biologi¬ 
cal and industrial processesi*^. Examples of polydis¬ 
perse systems are: plastic materials, liquid foams, as¬ 
phalt, micelles, and more. The property that typically 
varies among the constituent molecules is the charge, 
mass, and size. The fluids’ dispersity naturally induces 
a rich phenomenology not encountered in their single¬ 
component counterparts. As a consequence, polydisperse 
fluids have been the focus of a number of studies both 
theoretically^^— and via computer simulationsiir— , as 
well as in experiments^^Ti^. 

These investigations revealed that polydisperse flu¬ 
ids exhibit a very rich phase diagram with shifts in 
critical point and freezing temperatures, as well as 
many-solid/fluid coexistence regionsiML. The coexist¬ 
ing phases in fact show fractionation, i.e., the dispersed 
variable displays a distribution different from the par¬ 
ent in each phase. Moreover, polydisperse fluids ex¬ 
hibit terminal polydispersity beyond which the solid is 
not stabl o^°i^^i^^ . Polydispersity is thus regarded as an 
important factor to improve the glass-forming ability of 
liquids and to control the liquid fragilit y^^i^^ i^i^~— . 

Predicting the behavior of polydisperse fluids is a 
complex task, and several theoretical approaches have 
emerged to face the challenges of polydisperse fluids (see, 
e.g.. Refs, i i [TTI . and [2^ . The complex nature of 
the problem is enhanced by the fact that quasi-universal 
relations^^— for monatomic fluids, possessing great pre¬ 
dictive power, break down in increasingly polydisperse 
fluid a^^d^ . One such example is Rosenfeld’s excess en¬ 
tropy scalin g^^d^ in which a dimensionless transport co¬ 
efficient is correlated to the excess entropy with respect to 
an ideal gas. The latter observation for polydisperse flu¬ 


ids stands in contrast to confined fluids where the quasi- 
universal relations are very robust as the degree of con¬ 
finement is increased^^—. 

Explanations for quasi-universal relations 

in monatomic fluids were recently presented 
from a constant-potential-energy TV E17- dynamics 
perspective^^—, as well as from an exponential pair 
potential perspective^!^. In these perspectives, the 
quasi-universal relations are attributed to the monatomic 
fluids being so-called Roskilde-simple^^— . 

Roskilde-simple (RS) liquids have been detailed in a 
number of papers^^— and appear to be simpler than 
other types of liquids^. RS liquids are characterized 
by having strong correlations between virial and poten¬ 
tial energy equilibrium fluctuations at constant volume. 
In simulations of model liquids, the single-component 
Lennard-Jones (SCLJ) liquid, among many others, has 
been identified to belong to the class of RS liquids. Some 
experimental liquids have also been identihed as RS^“— . 
More specifically, the class of RS liquids is believed to 
include most or all van der Waals and metallic liquids, 
but to exclude most or all covalent-bonding, hydrogen¬ 
bonding, strongly ionic, and dipolar liquids^. Moreover, 
RS liquids are characterized by having isomorphs. Iso¬ 
morphs are curves in the phase diagram of RS liquids 
along which structure, dynamics, and some thermody¬ 
namic quantities are invariant in dimensionless units; a 
property that is not limited to the stable liquid region 
but also valid in the supercooled liquid and even out-of- 
equilibriu m^^d^ . 

Recently, Bagchi and coworkers^^*^ studied Linde- 
mann’s melting and Hansen-Verlet’s freezing criterion in 
size polydisperse LJ liquids. The authors showed sur¬ 
prisingly that terminal polydispersity can be analyzed 
in terms of these rules. It was also shown that order- 
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parameter maps for these liquids are not independent but 
highly correlated. We note that melting and freezing lines 
of RS liquids are isomorphs from which the rules of Lin- 
demann and Hansen-Verlet follow^. In fact, correlated 
order-parameter maps also follow from the existence of 
isomorphs in the phase diagrant^. 

The observations of Bagchi and coworkers imply that 
polydisperse LJ liquids retain the character of RS liq¬ 
uids. Conversely, the breakdown of quasi-universal rela¬ 
tions in polydisperse fluids may imply that polydisperse 
LJ liquids do not remain RS. Thus, we investigate in the 
following the effect of size polydispersity on strong virial- 
potential energy correlation in a systematic way to shed 
light on melting and freezing rules and quasi-universal 
behavior in polydisperse fluids. We note that previous in¬ 
vestigations of multi-component RS liquids focused only 
on binary mixture o^°i^^'^^ . 

Section|n]describes the applied methods and simulated 
model systems. Section Ell introduces RS liquids and 
explain their characteristics. In Sec. 113 we show re¬ 
sults from simulating size polydisperse LJ liquids in both 
two and three dimensions, in particular, the effect of size 
polydispersity on strong virial-potential energy correla¬ 
tion. Section 13 discusses the physics behind our observa¬ 
tions. Finally, Sec. I VII summarizes our work and presents 
an outlook. 


II. SIMULATION AND MODEL DETAILS 


We apply NVT molecular dynamics computer 
simulations2ii“2^ to study size polydisperse LJ liquids in 
two and three dimensions. The RUMD package (see 
http://rumd.org) is used for three-dimensional simula¬ 
tions while an in-house code is used for two-dimensional 
simulations. Pros and cons against applying constant- 
particle-number simulations to study polydisperse fluids 
are discussed in Ref. [zl. 

The polydisperse LJ liquid provides a convenient 
framework for initiating theoretical investigations of 
polydispersityiSi^, but has also been applied to mimic 
the behavior of for instance hydrocarbon or polymer so¬ 
lutions, idealized models for spherical micelles in water 
above the critical micelle concentration, and polystyrene 
lattices stabilized with polyoxyethylene chains^. 

If the total potential energy of the system is given by 
U = the size polydisperse LJ pair potential 

is 


_ 


Ary/ 

\rij) . 


( 1 ) 


where (Ty and e set, respectively, the length and energy 
scale of the pair interaction between particle i and par¬ 
ticle j {i,j = I, ...,iV, in which N is the number of par¬ 
ticles). As the study is focused on size polydisperse LJ 
liquids: e is independent of the pair interaction and fixed 
in simulations. The parameter (jy is defined from the 


Lorentz-Berthelot mixing rule2^ given by 


cTij = {ai + aj)/2. ( 2 ) 

In this study, ai and aj are evenly distributed on an 
interval from cr™™ to The polydispersity of the 

liquid S is defined by 


S = [pt) I = 


2 g.max _ ^min 

^max _j_ Q-min ’ 


(3) 


where and di are, respectively, the dispersed 

variable’s standard deviation and mean. We study poly- 
dispersities in the range 5 = 0% to J = 52% and apply a 
truncated-and-shifted pair potential cutoff at Tc = 2.5(Ty. 
We note that care should be taken when comparing 
the absolute polydispersity number 5 amongst different 
distributions!^. 


The units of the simulation are defined by setting the 
mean values over the N particles to unity, i.e., = 1, 

^NN — 1: = 1- -^11 particle masses are identical 

(and equal to unity). Figure [T] shows pair potentials 
of Eq. dD) with (Ty = 0.1, 1.0, 1.9 and e = 1. The 
latter corresponds to the range of LJ potentials used 
in simulations of a size polydisperse LJ liquid with 5 = 
52% (i.e., CT™" = 0.1 and = 1.9). For this high 
polydispersity the simulated pair potentials cover a very 
broad range of potentials. 



FIG. 1: The range of LJ potentials used in simulations of a 
polydisperse LJ liquid with <5 = 52%. 


The two-dimensional simulations use N = 1000 parti¬ 
cles, whereas the three-dimensional simulations use N = 
5000 particles. The finite-size effects are larger in two di¬ 
mensions than in three dimensions. In either case, how¬ 
ever, finite-size effects are observed to be small for the 
quantities studied in this paper and no additional aver¬ 
aging is performed. 
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III. ROSKILDE-SIMPLE LIQUIDS 


RS liquids and their associated simple properties have 
been detailed in a number of papers, starting with a series 
of five papersSSr—. RS liquids are characterized by hav¬ 
ing strong correlations between virial W and potential 
energy U equilibrium fluctuations at constant volume. 
The virial-potential energy correlation is quantified via 
Pearson’s correlation coefficient i? given by 


(AWAU) 

^{{AWY)./mJfy 


(4) 


in which A is deviation from mean values, (...) denotes 
NVT ensemble averages, and — 1 < i? < 1. RS liquids are 
defined pragmatically by requiring R > 0.90. R depends 
on the state point investigated but has been shown to be 
high in large parts of the phase diagram or not at all. 

Only inverse power-law (IPL) fluids with r“" pair po¬ 
tentials are perfectly correlating (i? = 1 since AW = 
{n/3)AU), but many models^i^i^^, such as the SCLJ 
liquid, the Kob-Andersen binary LJ mixture, the Lewis- 
Wahnstrom OTP model, the asymmetric dumbbell 
model, identical chain length LJ-bead polymer models, 
and some experimental liquids^^— , have been shown to 
belong to the class of RS liquids. 

The origin behind strong correlation in the SCLJ liq¬ 
uid is as follows. The LJ pair potential can be approx¬ 
imated by an IPL potential in the repulsive part^Zr— , 
i.e., below the potential mininum with an IPL-exponent 
n ~ 18. The required exponent is larger than the repul¬ 
sive term encountered in the LJ potential as the 
effect of attraction is to increase the steepness of the 
repulsive part. However, contributions to the fluctua¬ 
tions come from pair distances around the pair potential 
minimum^i^. As a consequence, one must add a lin¬ 
ear term and a constant to the IPL potential to obtain a 
proper description of the LJ potential: 


VLj{r) ^ Ar-^ + B + Cr. (5) 

The constraint of constant volume in the NVT ensem¬ 
ble has the following effect: When one nearest-neighbor 
distance increases another one decreases. Thus, upon 
summation the contribution from the linear term to U 
and W is almost constant and vanishes with respect to 
the fluctuations. For a very recent review of RS liquids, 
see Ref. HO- 

As mentioned, RS liquids are associated with a number 
of simple properties. RS liquids are, for instance, charac¬ 
terized by having isomorphs to a good approximation^. 
Consider two state points in a liquid’s phase diagram 
with density and temperature (pi, Ti) and (p 2 , T 2 ). 
These two state points are defined to be isomorphic if the 
following holds: Whenever a microconfiguration of state 
point (I) and of state point (2) have the same reduced 
coordinates, (R denotes a 3N- 


dimensional configurational-space vector and p = N/V, 
in which V is the volume), these two configurations have 
proportional Boltzmann factors, i.e., 

g-(7(RW)/feBTi ^ g-(7(R'"’)/feBT2^ (-0^ 

Here C 12 is a constant and depends only on the state 
points (I) and (2). An isomorph is defined as a continu¬ 
ous curve of state points that are all pairwise isomorphic. 
For IPL fluids C 12 = 1, and isomorphs thus encapsu¬ 
late the well-known scaling properties^ of these systems. 
However, for all other RS liquids C 12 ^ 1; accordingly, 
the state points can be represented effectively on a one¬ 
dimensional phase diagram only for certain properties^. 

An immediate consequence of isomorphs is that a num¬ 
ber of quantities are invariant along these curves. The 
invariants include the structure and dynamics in reduced 
units and some thermodynamic quantities, e.g., the ex¬ 
cess entropy with respect to an ideal gas Sex = Sex/N. 
We note that reduced units refer to macroscopic quanti¬ 
ties, such as and T, rather than the usual ap pro ach 

in simulations via microscopic quantities (see Ref. l63| for 
details). 

For RS liquids, temperature separates^ in a function 
of density and of excess entropy: ksT = h{p)f{sex)- 
From this equation, one observes the isomorph invariant 
h{p)/T which was used to define the so-called isomorph 
scaling^i^, in which the reduced relaxation time Tq, is a 
function fa = f{h{p)/T). Isomorph scaling explains clas¬ 
sical density scaling^^ for which fa = f{p'^/T), where 7 is 
a htting exponent. Density scaling is valid only for rela¬ 
tively small density variations, whereas isomorph scaling 
works even when the density variation is large^. Addi¬ 
tionally, the isomorph invariant h(p)/T provides a con¬ 
venient way to generate isomorphs since h{p) inherits the 
analytical structure of the potential energy function^. 

There exist additional simple observations for RS liq¬ 
uids than mentioned so far, e.g., the dominance of first- 
coordination shell interactions^!^, extension to non¬ 
trivial boundary conditions such as nanoconfinementZ^ or 
shear-flowZS, crystalline isomorphs^, and the relevance 
of Rosenfeld-Tarazona’s expressions^^ for these liquids^!. 
The interested reader is referred to these and other pa¬ 
pers for more information. 

Finally, a reformulation of isomorphs was very re¬ 
cently presented^. If Ra and Rf, denote two different 
3A-dimensional configurational-space vectors of a cho¬ 
sen state point, the isomorphic condition is 

U{Ka) < C/(Rfc) ^ U{XRa) < U{\Rb), (7) 

where A is a scalar. This definition encapsulates the pre¬ 
vious definition of isomorphs as a (very good) first-order 
approximation^. The mentioned invariants are exact 
within this reformulation, too. We note, however, that 
the constant volume heat capacity is not formally invari¬ 
ant. 
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IV. RESULTS FROM POLYDISPERSE LJ 
LIQUIDS 


The simulations in this paper were performed both in 
two and three dimensions, and thus we discuss each case 
separately beginning with the two-dimensional results. 



A. Simulations in two dimensions 


The investigation is initiated by considering the effect 
of size polydispersity on the structure and dynamics of 
the SCLJ liquid in two dimensions at p = 0.75 and T = 
0.70. For 0% polydispersity, this state point corresponds 
to the liquid phase. 

Figure Ufa) presents radial distribution functions 
(RDFs) upon the increase of polydispersity from S = 0% 
to ^ = 52%, whereas Fig. [IJb) shows bond-orientational- 
order function a^^'^^ (BOFs) defined via 



r 


^ 2^rAriV(JV - 1) (8) 

The division with g(r) in Fig. [2][b) approximately ac¬ 
counts for effects due to positional order in the liqui d^^i^° . 
In the equation above, <5 denotes Dirac’s i5-function, = 
Vi - rfc, L is the box length, and ipl = l/rn J2m=i 
in which is the number of nearest neighbors of par¬ 
ticle i, and 0^ is the angle between the vector r^i 
and the x-axis (particle m is a neighbor of particle i). 
Nearest-neighbor particles are identified via a Voronoi- 
construction^^. Finally, Fig. He) shows mean-square 
displacements (MSDs). 

Both measures of structural order (Figs. 

(b)) are seen to vanish as polydispersity is increased 
from 6 = 0% to i5 = 52% at constant density and 
temperature. The MSD shows approximately one 
order-of-magnitude decrease in the diffusion coefficient 
over the same variation in polydispersity. Introducing 
size polydispersity has thus a pronounced effect on 
the SCLJ liquid. The latter results confirm several 
previous simulations of polydisperse hard-sphere and LJ 
liquidsi^d^i^. We note that if the volume (area) fraction 
(j) = 1/L^ ’’■(o'j/2)^ of the system is kept constant, 

instead of density, a similar diminishment of structural 
correlations is observed. However, in the case of the 
dynamics an enhancement of the mobility is seer>2^. 



FIG. 2: Structural and dynamical effects of size polydispersity 
on the SCLJ liquid in two dimensions at p = 0.75 and T = 
0.70. (a) Radial distribution functions (RDFs). (b) Bond- 
orientational-order functions (BOFs, see Eq. (O). (c) Mean- 
square displacements (MSDs). Size polydispersity is seen to 
diminish structural order and slow down dynamics at constant 
density and temperature. 

Furthermore, Fig. |3] shows a snapshot of a particle 
configuration for the size polydisperse LJ liquid with p 
= 0.75, T = 0.70, 5 = 52%. Each particle is represented 
according to its value of Ui. A multitude of different 
sizes are present at this high polydispersity exemplifying 
again the rich and non-trivial diversity of the liquids 
studied. 






























FIG. 3: A snapshot of a particle configuration {N = 1000) for 
a highly size polydisperse LJ liquid with p = 0.75, T = 0.70, 
S = 52%. Each particle is represented according to its value 
of ai'. As the LJ potential is continuous some of the disks 
overlap slightly. A multitude of different sizes are present at 
this high polydispersity. 

The study is now turned towards the assessment of 
the effect of size polydispersity from the viewpoint of 
RS liquids. In this connection, we note that previous 
simulations of RS liquids were performed only in three 
dimensions. The simulations in this section are two- 
dimensional. 

FigureHlshows the correlation coefficient R as function 
of polydispersity S at both constant density p (black 
data points) and constant volume fraction cj) (red data 
points). For the three-dimensional SCLJ liquid, a 
typical value for the correlation coefficient is i? « 0.95. 
The inset shows the density-scaling exponen t^^i^^ 7 = 
{AUAW)/{{AU)'^) as function of polydispersity (see 
Sec. Imli . 

At constant density, we observe both a decrease and 
an increase in R depending on the starting state point at 
0% polydispersity. This observation is expected as the 
phase diagram for polydisperse liquids shows very rich 
phenomenologyiSii^. Although a decrease in R is ob¬ 
served for the two highest densities {p = 0.95 and p = 
0.85) all state points investigated are RS. For the lowest 
density {p = 0.75) the liquid becomes RS at high poly- 
dispersities. In fact, the results show that polydisperse 
LJ liquids at 52% polydispersity are RS liquids. This 
observation is non-trivial given the complex nature of 
polydisperse liquids (see also Fig. El). 

In the case of constant volume fraction, however, we 
observe a significant decrease in R as function of poly¬ 
dispersity. An increase in polydispersity, at constant vol¬ 
ume fraction, is accompanied by a decrease in pressure. 
This is in contrast to constant density which increases the 
pressure; Table |T] shows selected polydispersities and cor- 
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FIG. 4: Effect of polydispersity <5 on the correlation coefficient 
R (main figure, see Eq. 0 ) and density-scaling exponent 7 
= {AUAW)/{{AU)'^) (inset) for the polydisperse LJ liquid 
in two dimensions. Black data points are for constant density 
p, and red data points are for constant volume (area) fraction 
(j). For 0% polydispersity, p = 0.85 is inside the coexistence 
region, and p = 0.95 is inside the crystalline region. 

responding volume fractions along the three simulated 
constant-density curves. For the density-scaling expo¬ 
nent 7 , we only observe a decrease in its value, when 
polydispersity is increased. 


5 

<). (p = 0.75) 

0 (p = 0.85) 

(p = 0.95) 

0 

0.589 

0.668 

0.746 

11.55 

0.597 

0.676 

0.755 

23.09 

0.619 

0.702 

0.784 

34.64 

- 

0.746 

0.834 

51.96 

0.746 

0.845 

- 


TABLE I: Selected polydispersities J and corresponding vol¬ 
ume (area) fractions (fi along the three simulated constant- 
density curves. 

As mentioned previously, the quantity h{p)/T is in¬ 
variant along an isomorph (see Sec. IIII|) . The function 
h(p) for a LJ liquid in two dimensions is given by 

M/5) = (7*/3-1)/5® + (2-7*/3)p^ (9) 

In the equation above: p = p/p*, p* is a chosen reference 
density, and 7 * is the value of 7 = (AC/AFF)/((AC/)^) ob¬ 
tained from the equilibrium fluctuations at p* (see Refs. 

and m for additional details). In this and the follow¬ 
ing section the reference state point is p* = 1 and T, = 
1. We now proceed to study isomorphs of a moderately 
polydisperse LJ liquid. Isomorphs are generated by the 
procedure stated belo w^^i^^ : 

1. A starting state point is chosen. 

2. Density is varied by some percentage (typically 10- 
20 % but can be even a factor of two^). 
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3. The temperature of the isomorphic state point is 
calculated via Eq. m keeping h{p)/T constant. 

4. A simulation is performed at the predicted state 
point and the procedure is repeated. 

Figure EJa) shows RDFs along an isomorph with 20% 
density increase ( 7 * = 7.02 and starting state point 
p = 0.85 and T = 0.70). For reference, Fig. E^b) 
shows an isotherm with the same density increase. A 
good invariance of structure is seen along the isomorph, 
whereas structure on the isotherm clearly shows a poorer 
scaling. 





FIG. 5: RDFs for a moderately polydisperse LJ liquid (<5 
= 28.87%) in two dimensions, (a) Along an isomorph. (b) 
Along an isotherm. The highest density state point is not 
fully equilibrated. 

Figure El shows the corresponding BOF figures. A 
good invariance is again seen along the isomorph but 
not along the isotherm. 

Turning to the dynamics in Fig. (3 a near-perfect 
collapse is seen for the MSD on the isomorph whereas, as 
expected, a poor invariance is seen along the isotherm. 
The fact that liquids are highly size polydisperse is thus 
not a hindrance to showing simple isomorphic behavior. 



FIG. 6: BOFs for a moderately polydisperse LJ liquid {5 
= 28.87%) in two dimensions, (a) Along an isomorph. (b) 
Along an isotherm. The highest density state point is not 
fully equilibrated. 


B. Simulations in three dimensions 


This section studies size polydisperse LJ liquids 
in three dimensions. Figure [ 8 ] shows the correlation 
coefficient R as function of polydispersity at constant 
density and volume fraction. As for the two-dimensional 
results both a decrease and an increase in R is observed, 
when polydispersity is increased at constant density 
(i.e., increasing pressure). For the density p = 0.75, the 
liquid becomes RS at high polydispersities, whereas for p 
= 0.85 and p = 0.97 all state points investigated are RS, 
albeit R decreases. We observe again that highly size 
polydisperse LJ liquids are RS: a non-trivial observation. 

























































FIG. 8 : Effect of polydispersity <5 on the correlation coefficient 
R (main figure, see Eq. ®) and density-scaling exponent 7 = 
(AUAW)/{{AU)'^) (inset) for the polydisperse LJ liqnid in 
three dimensions. Black data points are for constant density 
p, and red data points are for constant volume fraction cj>. 
The highest polydispersity data points of p = 0.75 and p = 
0.85, and the two highest polydispersity data points of p = 
0.97 are not fully equilibrated due to very slow dynamics. For 
0% polydispersity: p = 0.97 is inside the crystalline region. 
The magenta asteriks data points span densities (T = 0.70): 
0.97, 0.85, 0.75, 0.70, 0.65, and 0.60. 



FIG. 7: MSDs for a moderately polydisperse LJ liquid (J 
= 28.87%) in two dimensions, (a) Along an isomorph. (b) 
Along an isotherm. The highest density state point is not 
fully equilibrated. 


In the case of constant volume fraction, a significant 
decrease in R is observed and is accompanied by the virial 
turning negative (but not the pressure). Actually, the 
liquid goes from RS to non-RS in the studied range of 
polydispersities. For single-component liquids: R usually 
decreases rapidly, when the virial becomes negative^. 

The density-scaling exponent 7 shows — as in two di¬ 
mensions — a decrease in its value as function of poly¬ 
dispersity. Actually, it is observed from Fig. [5] that all 
constant-density curves collapse at high polydispersities 
with 7 Ri 4.5. Even a path at 5 = 52%, where density 
varies from p = 0.97 to 0.60 {(j) = 0.94 to 0.58) and the 
virial from W = 110 to 5.8, gives approximately the same 
value of 7 (see magenta asteriks data points). 

This observation indicates that for highly size polydis¬ 
perse LJ liquids 7 is mainly controlled by the degree of 
polydispersity S. In the non-reformulated isomorph the¬ 
ory 7 is predicted to be a function only of density^ (see 
Sec. mil) . For the SCLJ liquid, a value of 7 ~ 5.5 is ob¬ 
tained at low-to-moderate pressures (see inset and Ref. 
[^. 

Very recently, an expression was established relating 
7 to the pair potential of monatomic RS liquids^!. The 
expression is given by 


7(p) = 
n^^^(r) = 


(r) 



v(^)(r) ’ 


( 10 ) 

( 11 ) 


where A is a dimensionless number, and is the p’th 





















derivative of the pair potential f. 7 is then related to the 
curvature and torsion of the pair potential. A is chosen 
as the reduced distance (f = rp^/^) at which r^gir) ob¬ 
tains its maximum, i.e., the most likely nearest-neighbor 
distance^i. In this case, however, 7 is not a function only 
of density. 

RDFs along an isoinorph and isotherm for a polydis- 
perse LJ liquid with (5 = 23.09% are given in Fig. |9] ( 7 * 
= 4.79 and starting state point p = 0.97 and T = 0.70). 
The isomorph is generated via the aforementioned proce¬ 
dure keeping h{p)/T constant. The expression^ for h{p) 
in three dimensions is given by 


( 13 ) 

in which 


h(p) = (7*/2-l)p4 + (2-7,/2)p2, (12) 

with definitions provided in the previous section. The 
RDF is to a good approximation invariant along the 
isomorph and less so on the isotherm (the highest 
density state point is, however, not fully equilibrated 
due to slow dynamics). 




The above expression for gives coarse-grained 

Steinhardt bond-orientational-order parameter o^°'^^ . 
Igm is the spherical harmonic function with degree I = 
6 and order m = [— 6 , 6 ]. The notation rii + 1 indicates 
that the sum is taken over all the nearest neighbors of 
particle i including particle i itself; nearest-neighbor 
particles are identified via the twelve-closest particles. 
From Fig. uni a similar conclusion as for the RDF 
is reached; a good invariance is observed along the 
isomorph and less so on the isotherm. These results 
also indicate a good invariance of the first coordination 
shell (FCS) environment along the isomorph and are 
consistent with the role of the FCS in determining the 
structure and dynamics of RS liquids (see Ref. [63 for 
more information on FCS results). 



FIG. 9: RDFs for a moderately polydisperse LJ liquid {S = 
23.09%) in three dimensions, (a) Along an isomorph. (b) 
Along an isotherm. The highest density state point is not 
fully equilibrated. 

Figure uni shows BOFs along the same isomorph and 
isotherm. The BOF is defined in three dimensions via 
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r (reduced units) 

FIG. 10: BOFs for a moderately polydisperse LJ liquid (<5 
= 23.09%) in three dimensions, (a) Along an isomorph. (b) 
Along an isotherm. The highest density state point is not 
fully equilibrated. 


Finally, the dynamics in terms of the MSD is shown 
in Fig. [TT] An excellent invariance is obtained along the 
isomorph, whereas the isotherm shows approximately 
five orders-of-magnitude reduction in the diffusion 
coefficient. 




FIG. 11: MSDs for a moderately polydisperse LJ liquid (S 
= 23.09%) in three dimensions, (a) Along an isomorph. (b) 
Along an isotherm. The highest density state point is not 
fully equilibrated. 


V. DISCUSSION 

We have demonstrated that both two- and three- 
dimensional size polydisperse LJ liquids belong to the 
class of RS liquids up to a significant level of polydisper- 
sity. This shows that the theory of isomorphs can be ap¬ 
plied to multi-component systems which contributes to 
the simpler physical description of this important class 
of liquids. Nevertheless, our simulations also showed un¬ 
usual behavior when compared to the SCLJ liquid. 1) 
At higher densities, E decreased with increasing poly- 
dispersity (pressure). An increase in pressure normally 
reduces interparticle distances in which a pure IPL de¬ 
scription is a better approximation, and one would a 
priori expect R to increase. 2) For high polydispersi- 
ties, the density-scaling exponent 7 is controlled only by 
the degree of polydispersity exhibiting also significantly 
lower 7 -values at low-to-moderate pressures (7 4.5 ver¬ 

sus JSCLJ ~ 5.5 in three dimensions). Below we try to 
understand these observations in more detail presenting 
three-dimensional results only. 

Figure [12] shows like-like RDFs (colored curves) for the 
small, medium, and large particles at p = 0.85, T = 0.70, 
S = 52%, and the corresponding SCLJ liquid RDF (black 
curve). The RDFs for the medium and large particles 
move to shorter distances when compared to the SCLJ 
liquid (note the scaled distance axis). Conversely, the 
RDF for the small particles moves to larger distances 
and shows a very broad distribution. The higher value 
of the first peak is due to a clustering tendency of small 
particles in spaces surrounded by bigger particles (see 
Fig. |3|) 

The results signify a decoupling, with respect to pres¬ 
sure, amongst the small and large particles and could be 
the reason why R decreases with increasing polydisper¬ 
sity (pressure) at higher densities. To verify this, we sim¬ 
ulated at the same state point (not shown) a truncated 
uniform distribution excluding all the small(er) particles 
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(i.e., particles with cr^ < 1). In this case: R increased 
with polydispersity (pressure). 

Despite that the small particles are located around 
the minimum of the LJ potential, fluctuations in their 
positions are quite large. We speculate that this large- 
amplitude positional fluctuation might diminish the 
impact of the particle-size dependent inhomogeneous 
stress distribution on R (see Eq. (O). Accordingly, it 
may also explain why the liquid remains RS up to such 
high polydispersities. Additional studies are, however, 
needed to clarify this issue. 



FIG. 12: Like-like RDFs for the small (s), medium (m), and 
large (/) particles at p = 0.85, T — 0.70, <5 = 52% (colored 
curves; W = 55.4) in three dimensions. gss(r) uses particle 
sizes in the range (Ji = [0.1, 0.2], gmm{r) uses particle sizes in 
the range Oi = [1.0, 1.1], and gu{r) uses particle sizes in the 
range at = [1.8,1.9]. We note that this state point is not fully 
equilibrated. The black curve gives SCLJ liquid RDF at p = 
0.85 and T = 0.70 {W = 0.31). The magenta dashed curve 
gives the LJ pair potential, and the black vertical dashed line 
delimits the minimum of the LJ potential. 


The investigation is now focused on 7 -values at low-to- 
moderate pressures for high polydispersities (7 ~ 4.5). In 
the light of the theoretical results on monatomic RS liq¬ 
uids (see Eq. (ITUI) !. Fig. fOI al shows r^gaa{r) for the 
small (s), medium (m), and large (/) particles {a = s, m, 
0 at p = 0.60, T = 0.70, 5 = 52% (W = 5.8). Figure 
[T5f bl shows rfi^r) for the LJ potential (Eq. (ITT]) 1. Dis¬ 
regarding now the fact that the liquid is polydisperse we 
apply (r) to make predictions for 7 ^ using the proce¬ 
dure mentioned in Sec. lIVBI for RS monatomic liquids. 

In this way 7 ^ = 5.5 and 7 ; = 4.9 are obtained for 
the medium and large particles, respectively, whereas 
7 s = —2.9 gives a negative value due to the divergence in 
Interestingly, 7 ^ gives a value very close to that 
expected for a SCLJ liquid at low-to-moderate pressures 
ilSCLj ~ 5.5). Nevertheless, the results indicate that 
the large particles play an important role in determining 
the value of 7 for polydisperse LJ liquids with 7 ; = 4.9. 
It is, however, unclear which role the small particles 
play, as r^gssir) peaks at very large distances with a 
resulting negative 7 ^. We speculate that there is some 


cancellation between the medium particles of large 7 
and the small particles of negative 7 . More in-depth 
investigations are, however, needed. 



r/o.. 

ij 


(b) 


9 . 0 - 
8.0 - 
7.0 - 
6.0 - 
5 . 0 - 
4.0 - 


Equation {11) for the LJ potential 


I _ 


0.6 
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FIG. 13: (a) r^gaa(r) for the small (s), medium (m), and 
large (1) particles (a = s, m, 1) at p = 0.60, T = 0.70, J = 
52% {W = 5.8) in three dimensions. gss{r) uses particle sizes 
in the range at = [0.1, 0.2], gmm{r) uses particle sizes in the 
range Ui = [1.0, 1.1], and gii(r) uses particle sizes in the range 
(Ti = [1.8, 1.9[. (b) 'nS^\r) of Eq. dill) for the LJ potential. 


VI. CONCLUSION AND OUTLOOK 

Polydisperse fluids show rich phenomenology with re¬ 
spect to their single-component fluid counterparts and 
manifest themselves in as diverse systems as bitumen in 
road construction^^ to micelles in biological systems^. 
In spite of this fact, we have shown that even highly size 
polydisperse LJ liquids belong to the class of RS liquids. 
As a consequence, polydisperse LJ liquids inherit the as¬ 
sociated simple properties of RS liquids^i^ and are ex¬ 
emplified here by studying isomorphs of moderately poly¬ 
disperse LJ liquids in both two and three dimensions. 

The SCLJ liquid is a prime example of a RS liquid, 
however, many more model and experimental liquids are 
identified to belong to this class of liquids. A natural 
question thus arises: How are the remaining RS liquids 
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affected by introducing size polydispersity? As an ex¬ 
ample, identical chain length LJ-bead polymer models 
have very recently been identified as RS^. For polymers 
the relevant size polydispersity is given by a distribu¬ 
tion of chain lengths (mass distribution); can one expect 
such polymeric liquids to be RS? The answer is most 
likely in the affirmative. It is known from experiments 
that scaling relations typical of RS liquids, such as den¬ 
sity scaling and excess entropy scaling, are also valid for 
polymers^. As another example, rigid-bond molecular 
liquids with LJ interactions have also been identified as 
RS^; it is expected that these liquids remain RS when 
introducing size polydispersity. A more general hypoth¬ 
esis is that all RS liquids will remain so to a high degree 
of accuracy when introducing size polydispersity. Future 
research should, however, focus on clarifying to what ex¬ 
tent this hypothesis is valid and on understanding the 
implications thereof. 

To conclude, the theory of isomorphs extends readily 
to multi-component systems, and when doing so it can 
be applied to improve even further the understanding of 
these intriguing systems. Experimental liquids are of¬ 
ten distributed in additional or other variables such as 
mass or charg o^d^d^ . More work in this direction is also 
needed to understand how these factors influence the re¬ 
sults presented here. 
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